Spin-torque shot noise in magnetic tunnel junctions 
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Spin polarized current may transfer angular momentum to a ferromagnet, resulting in a spin- 
torque phenomenon. At the same time the shot noise, associated with the current, leads to a 
non-equilibrium stochastic force acting on the ferromagnet. We derive stochastic version of Landau- 
Lifshitz-Gilbert equation for a magnetization of a "free" ferromagnetic layer in contact with a "fixed" 
ferromagnet. We solve the corresponding Fokker-Planck equation and show that the non-equilibrium 
noise yields to a non-monotonous dependence of the precession spectrum linewidth on the current. 



Magnetization dynamics of a ferromagnet under influ- 
ence of a spin polarized current is a subject of intensive 
investigations (for recent reviews see Rcfs. [2,0). It 
was realized 0, Q that the spin current may transfer the 
angular momentum to the ferromagnet, resulting in a 
torque acting on its magnetization direction. In the case 
of a small ferromagnetic domain the torque may lead to a 
rotation of the magnetization as a whole, rather than to 
an excitation of spin waves. This phenomenon, allowing 
for an electronic manipulation of the magnetization, has 
a promise for a number of potential applications. 

The effect has been recently observed 0, H, 0, H m a 
setup, where the spin-torque magnitude and direction are 
tuned to compensate exactly the dissipation force acting 
on the magnetization of the "free" ferromagnetic layer. 
This leads to an undamped precession which is detected 
through the induced microwave radiation. Both the spec- 
tral width and the generated power exhibit a strong de- 
pendence on the current flowing through the interface of 
the two ferromagnets. It was shown later that the 

equilibrium thermal noise, first introduced in dynamics of 
micromagnets by F-L. Brown [lij ]. may partially account 
for the observed linewidth. 

On the other hand, since the experiments are per- 
formed under non-equilibrium conditions (spin current 
strong enough to balance the dissipation), one needs to 
address non-thermal sources of noise as well. The most 
essential of them is the spin shot noise associated with the 
discreteness of spin passing through the interface. This 
phenomenon may be accounted for by adding a fluctuat- 
ing part to the spin current vector in the Slonczewski's 
torque term of Landau-Lifshitz-Gilbert (LLG) equation 
l s — * I s + SI s (t). The resulting stochastic LLG equation 
for the unit vector m = M/M in the direction of the 
magnetization M takes the form 



— = -7 m x H cff + a{6) 
at 



m x 



dm 
~~dt 



7 
MV 



[m x [(I s + Sl s ) x m] ] 



(1) 



Here 7 is gyromagnetic ratio, H c s = —dF/dM. is the 
effective magnetic field, which includes both an external 
field and magnetic anisotropy, and V is a volume of the 



free ferromagnet. Gilbert damping a(9) is renormalized 
by the coupling to the fixed ferromagnet 0, [l3| and is 
thus dependent on a relative orientation angle 8 of the 
fixed and free ferromagnets. 

One could expect that the fluctuating part of the spin 
current vector Sl s (t) is preferentially directed along the 
spin polarization of the incoming electron flux, i.e. along 
I s . This is not the case, however, due to the quantum 
nature of the effect. Indeed, each spin-flip event transfers 
exactly one h unit of the angular momentum to the free 
ferromagnet. Due to the uncertainty principle, direction 
of an ensuing magnetization rotation is completely ran- 
dom. As a result, the fluctuating part of the spin torque 
must have an isotropic correlator 



{Sl Sti Sl Stj ) = W{6)5 i:j 5{t-t'), 



(2) 



where the variance T>(0) does not depend on the cartesian 
indexes i, j — x,y,z, but may depend on the mutual ori- 
entation of the two ferromagnets. Because of the isotropy 
of the stochastic torque, it can be equally well repre- 
sented by a fluctuating magnetic field H fj — -> H c ff + h(t), 
instead of the fluctuating spin flux Sl s (t). In the lat- 
ter case the correlator of the stochastic fields reads as 
(hi(t)hj(t')) = 2V(MV)- 2 6ij5(t - t'). This type of non- 
equilibrium noise was considered in Ref. [l^ | in context 
of NFN structures. 

Here we consider a model of a magnetic tunnel junction 
(MTJ) Ref. 14|. Magnetization dynamics of the free 



ferromagnet is described using Holstein-Primakoff (HP) 
parametrization (lBj of its total spin operator by deriv- 
ing semiclassical equations of motions for HP bosons. We 
employ Keldysh formalism to allow for non-equilibrium 
conditions, i.e. voltage bias between the two ferromag- 
nets IitI . Il8j |. After integrating out the fermionic 
degrees of freedom in the second order in both tunnel- 
ing and spin-flip processes, we obtain an effective ac- 
tion for HP bosons, which encapsulates deterministic 
forces (external magnetic field and spin-torque) along 
with the stochastic term. The latter contains an informa- 
tion about both equilibrium and non-equilibrium noise 
components. 

The result of the program, outlined above, is the fol- 
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lowing noise correlator 



V(6) 



MV 

1 



■ a T ■ 



(3) 



where cto is a bare Gilbert damping of an isolated grain 
and V is a voltage bias between the two ferromagnets. 
The non-equilibrium part of the noise is proportional to 
the spin-flip current I s / (6) . The latter counts the number 
of real spin flip events (as opposed to the spin current I s 
which is associated with virtual spin flips). In the MTJ 
setup we found for the spin-flip conductance 



dv 



h 

4e 



(■//-> snr i - 



(4) 



Gp = Gj 



G — ; Gap — G-\ — + G. 



+ > 



where we adopted notations of Ref. for the par- 
tial conductances G aa > between the spin-polarized bands 
of the two ferromagnets. Here Gp/ap stay for elec- 
tric conductances in parallel (P) and antiparallel (AP) 
configurations, with Gp > Gap- The electric conduc- 
tance of the MTJ in an arbitrary orientation is given by 
dl e (0)/dV = G P cos 2 ((9/2) + Gap sin 2 (0/2). Notice that 
the spin shot noise is minimal for P orientation and max- 
imal for AP one - exactly opposite to the charge current 
and the charge shot noise. In the same notations the spin 
conductance is [3] 



dl s h 
dV ~ 4e 



(G H 



G__ 



G+- - G ) , 



(5) 



where the angular dependence is already explicitly taken 
into account in Eq. ([T]). A non-polarized current, i.e. 
G +(J = G-o-, does not exert deterministic spin torque, 
nevertheless it still induces the shot noise torque on the 
ferromagnet. This effect was recently discussed in the 
context of NFN structures (TqI ] . 

The same calculation also leads to a renormalization 
of Gilbert damping coefficient in LLG equation fTJ) 



a(9) = a 



tvy (dl s} {9) 



eMV 



dV 



(6) 



where the spin- flip conductance is given by Eq. ([J). 
Such an enhancement of the dissipation due to the spin 
transport between the two ferromagnets is discussed in 
Refs. @, EJ]. In compliance with the fluctuation- 
dissipation theorem, the equilibrium (V — > 0) noise cor- 
relator is given by V = TMVa(9)/j 0. Due to the 
renormalization, the damping is minimal in P orienta- 
tion. The angular dependence of the damping term se- 
lects a unique angle, where the spin-torque compensates 
dissipation. Upon increasing the spin current above the 
critical one, such an angle gradually increases from zero 

to 7T. 



Below we restrict ourselves to a setup, where both the 
magnetic field and the spin current are directed along 
the z-axis specified by the magnetization direction of the 
fixed ferromagnet. Transforming Eq. ([!]) to the spher- 
ical coordinates [2(|, one obtains two coupled Langevin 
equations for the relative angle 9(t) and the azimuthal 
angle cj)(t) 

9 = -sin0T 2 (0) + cot 9V{9) +&>(*), (7) 
</> = n(0)+csc0&(i); (8) 

T z (0) = i (a(6)H z + I 8 /(MV)) ; U = 1 H Z - aT z , 

where 7' = 7/(l+a 2 ) and V(9) = 2>(0hY / (MV) 2 . Here 
l;o(t) and £<j>(t) are two uncorrelated random noises with 
the correlators 



<&(*)&(*')> = <^(W)> = me)6(t - 1') 



(9) 



Both Eq. (TTJ) and Eqs. 0, © should be interpreted in 
the sense of retarded regularization, or Ito calculus [2(ij ]. 

Let us first analyze deterministic dynamics described 
by Eqs. (0, © with V -> 0. For a strong enough 
(and negative for positive H z ) spin current the condition 
T Z {S) = may be satisfied for a certain angle 8. Notice 
that it is the angular dependence of the enhanced Gilbert 
damping 0, , which is responsible for the angle selec- 
tivity. In such a case Eq. {5| describes a stable undamped 
precession with the frequency f2(0) = 7-Hz- The inten- 
sity of the induced microwave radiation 0, @, 0, S, 0] is 
given by the square of the oscillating magnetic moment, 
i.e. proportional to sin 2 8. 

To analyze effects of the noise we shall assume that 
a(8) -C I, allowing for time scales separation. The fast 
variable is the azimuthal angle <fi(t), while the angle 9 if) 
is the slow one. For a fixed slow variable 8 Eqs. JHJ), ([9|) 
lead to the Lorentzian shape of the emitted microwave 
power spectrum 



S(w, 8) oc sin 2 



2csc 2 6»2?(6») 



[lo - ni&)} 2 + csc 2 8T>(8) 



(10) 



To compare with the observed power spectrum this ex- 
pression should be averaged over the stationary prob- 
ability distribution P{9) of the slow degree of freedom 
S(w) = J sin8d8 P{8)S(u,8). The distribution function 
P(6,t) obeys the Fokker-Planck equation which follows 
from Eqs. ©, © [UlU 



P 



1 



sin ( 



sin 2 9T Z (9)P + sin 9 d e (v{9)P^j . (11) 



The stationary solution of Eq. (jlip is given by 



P{9) 



ZV(9) 



exp 



T ( 

sin 8' d8' z{ 



V(9') 



(12) 
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where constant Z is chosen to satisfy the normalization 
condition J Q sin 9 d9 P{&) = 1. For a weak noise the sta- 
tionary distribution function has a sharp maximum close 
to the angle 9, where the deterministic spin-torque com- 
pensates the dissipation. 

Figure Q] shows the calculated spectral linewidth at the 
half maximum as a function of the applied voltage. The 
origin of the non-monotonous dependence may be un- 
derstood by inspection of Eq. (fTU|) . The initial decline 
is due to the geometric factor esc 2 9 in the width of the 
Lorentzian. As the voltage increases, so does the angle 
where the distribution function exhibits the maximum. 
Due to esc 9, coming from the noise term on the r.h.s. of 
Eq. dHJ), the amplitude of the phase noise decreases lead- 
ing to a narrowing of the power spectrum. The initial de- 
crease of the spectral width was discussed in Ref. [l0|,[ll| 
in terms of the equilibrium thermal noise 

The subsequent increase of the spectral width at yet 
larger voltages was observed in a number of experiments 
0,13, S] and, to the best of our knowledge, remained 
unexplained. It naturally comes about due to the non- 
equilibrium component of the noise. Indeed the noise 
correlator ([3]) grows with the applied voltage due to the 
growth of the spin- flip current. The dependence of the 
spin- flip current on the bias is faster than the lin- 

ear because of the increase of the spin-flip conductance, 
Eq. 0, with the operation angle 9 on top of the over- 
all proportionality of L/ to the bias V. As a result, for 
eV > 2T the noise variance T> grows rapidly, leading to 
the broadening of the power spectrum, cf. Eq. (fT0|) . 

Another consequence of the non-equilibrium noise is 
the saturation of the spectral linewidth at small tem- 
peratures [9(. Since the devices are always operated at 
currents larger than the critical one, the noise intensity 
© is finite V{9) = (fi/2)I s/ (0) even at T -> 0. Thus de- 
creasing the temperature one should observe saturation 
of the linewidth at T ~ eV, provided the induced damp- 
ing (the second term on the r.h.s. of Eq. ©) is larger 
than the bare one. 

In the remainder of the paper we outline derivation of 
Eqs. © - ©■ The MTJ is modeled by the two itiner- 
ant ferromagnets, whose majority (a — +) and minority 
((j = — ) bands are described by the operators cj^, Cfc CT for 
the fixed ferromagnet and d\ , di a for the free layer. We 
found convenient to work in the instantaneous reference 
frame, where the magnetization of the free layer points 
in the z-direction. The corresponding Hamiltonian takes 
the following form 

Hq = '^2 tkaC^Ckc + ^(e; - JS z (j)d\ a di a - 7S H 



la 



-J (S+s_ + 5_s_| 



kl,aa ; 
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FIG. 1: (Color online) Calculated linewidth of the microwave 
spectral power vs. applied voltage in units of temperature. 
Inset: calculated microwave power (solid line) and the av- 
erage precession angle (dashed line) vs. voltage. Parame- 
ters: MV/{hi) = 10, hjH z = 3K, T = IK, a = 0.01. 
Conductances in units e 2 /h: Gp = 0.181, Gap = 0.019. 
dI s /dV = O.Ole. 



is (a er> 



<j\(t')W, where the spin-transformation matrix 
= e- i<T ^ 2 cos(9/2 and (a\a') = e 1 ' 



and S = I Y,laa> d 



and (ojcr') = e^/ 2 sin 0/2; 
u cto ./ uj CT ' is the spin of itinerant elec- 
trons, while s± = s x ± is y . We have explicitly accounted 
for the interactions of the itinerant electrons in the free 
layer with its total spin S = MV/7. To make the latter 
a dynamical variable we use HP parametrization [l5| 



S„ = S- b%; S- = bW2S-tfb; S+ = \/2S-tfbb, 

where b\ b are usual bosonic operators. 

Next we write the corresponding action in terms of 
complex fermionic and bosonic fields Ck<r(t), di a (t),b(t), 
where the time variable runs along the closed Keldysh 
contour 



17, m 



here the spin-dependent tunneling matrix elements are 



We then transform to two- 
component vector notations in terms of symmetric (clas- 
sical "cZ") and antisymmetric (quantum "q") combina- 
tions of the forward and backward propagating fields. 
One should keep in mind that the distribution functions 
of the c and d fermions have a relative shift of the chem- 
ical potentials by eV. The fermionic fields may be inte- 
grated out exactly and the remaining bosonic effective 
action expanded to the second order in the tunneling 
amplitude W and to the first and second orders in the 
spin- flip processes S±s^. The corresponding processes 
are represented by the diagrams of Fig. [2j The approx- 
imations are justified by the weakness of tunneling and 
largeness of S » h. 

The resulting action for the complex bosonic fields 
b c i{t), b q (t) takes the form S = Sq + S± + S2 where the 
subscript indicates the order in spin-flips processes. Here 
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FIG. 2: Diagrams for spin-flip processes: (a) the first order, 
describing the spin-torque; (b) the second order, describing 
the spin shot noise along with the enhanced damping. Solid 
(dashed) lines denote electronic propagators in the free (fixed) 
layers. Bold dashed lines are propagators of HP bosons (spin- 
flips). Tunneling vertices are denoted by circles with crosses. 



the bare action is [21 1 

S = Jdtb q {t) (id t b c i{t)+ 1 ^Sj2H+^ +c.c. (13) 

The first order correction in spin-flip amplitude is repre- 
sented by diagram of Fig. [2^. This is a virtual transition 
into an opposite spin band with a subsequent tunneling 
out of the free layer into the " correct" spin band of the 
fixed magnet. The latter process is possible for 9 ^ 0, n, 
due to a finite Wgf . The net result is transferring an- 
gular momentum h to the total spin of the free layer, i.e. 
the deterministic spin-torque [H, 0] • The corresponding 
contribution to the action is 



Si 



dt b q (t)I s sin Oe-^ + cc, (14) 



where I s is given by Eq. © with G aa , = ^-\Wfv a c 
and v° d are densities of states of the two ferromagnets 
in the a band. The second order processes in spin-flips 
are depicted by the diagram of Fig. [2)3. These are real 
(i.e. Golden rule) processes, which matrix elements in- 
clude the spin- flips. They lead to dissipation as well as 
fluctuations. The corresponding action is 



S-x = dt 



a(6)(b q d t b cl - b cl d t b q ) + ^V(6)b q b q 



,(15) 



field (51+ (t) =I S , X 
in b q (t) and b q (t). 



where T>{6) and a{6) are given by Eqs. ©, g]) and © 
(without internal dissipation ag)- 

One then decouples the last term on the r.h.s. of 
Eq. (fTS"]) by means of the complex Hubbard-Stratonovich 
h il s ,y- The remaining action is linear 
It constitutes thus resolution of func- 
tional ^-functions of the first order Langevin equations 
on b c i(t) = y/MV/(2'y) m + (t) and its complex conjugate. 
Those are nothing but m± components of Eq. |T]), with 
the noise intensity given by Eqs. ©-([I]) [21 1. 

To conclude: in the presence of the spin-torque, caused 
by a spin polarized current, the LLG equation should be 
modified to include a stochastic Langevin term which 



accounts for the shot-noise, associated with the spin cur- 
rent. This term is different from the previously dis- 
cussed thermal stochasticity in LLG equation [l2|, be- 
cause of its non-equilibrium origin. We have derived the 
corresponding noise correlator in the MTJ setup. We 
have argued that the non-equilibrium noise manifests it- 
self in a non-monotonous voltage dependence and low- 
temperature saturation of the linewidth of the precession 
power spectrum. 
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